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We study the occurrence of shell crossing in spherical weakly charged dust collapse in the presence of 
a non- vanishing cosmological constant. We find that shell crossing always occurs from generic time- 
symmetric regular initial data, near the center of the matter configuration. For non-time-symmetric 
initial data, the occurrence — or lack thereof — of shell crossing is determined by the initial velocity 
profile, for a given mass and charge distribution. Physically reasonable initial data inevitably leads 
to shell crossing (near the center) before the minimum bounce radius is reached. 

PACS numbers: 04.20.Dw, 04.40. Nr, 04.70.Bw 
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I. INTRODUCTION 

Tolman-Bondi metrics |l],|| describe the gravitational 
collapse of spatially bounded spherical dust configura- 
tions in an otherwise empty spacetime. The model con- 
sists of a general spherically symmetric metric, matched 
to a Schwarzschild exterior, with the matter content be- 
ing general inhomogeneous dust. Spherical symmetry im- 
plies that the initial data — two arbitrary functions giving 
the initial density and velocity profiles — are a function 
of the radial coordinate r alone. Since the collapse is 
pressureless, every shell of dust with finite proper radius 
will collapse through its Schwarzschild radius, become 
trapped and proceed to become singular at the center of 
symmetry in a finite amount of proper time (as measured 
by an observer comoving with the shell). 

For homogeneous dust distributions — Oppenheimer- 
Snyder collapse || — all the shells become singular at the 
same time, and thus none of the shells cross g. For in- 
homogeneous matter configurations, however, the proper 
time for collapse depends on the (comoving) coordinate 
radius r, and thus the piling up of neighboring matter 
shells at finite proper radius can occur, thereby produc- 
ing two-dimensional caustics where the energy density 
and some curvature components diverge [g|. These sin- 
gularities can be locally naked, but they are gravitation- 
ally weak [pLM — curvature invariants and tidal forces re- 
main finite — and, from the viewpoint of geodesic com- 
pleteness, analytic continuations of the metric can always 
be found, in a distributional sense, in the neighborhood 
of the singularity M. In this respect, shell crossings are 
not genuine physical singularities; rather, they signal the 
intersection of matter flow lines at a certain spacelike sur- 
face (in spherically symmetric geometries), and thus the 
breakdown of the model beyond that surface. (It is worth 
noting that they also occur in spherical inhomogeneous 
Newtonian gravitational collapse). 

The inclusion of an electric charge density in spheri- 
cal dust models has been considered by various authors 
P-|ll| ■ The principal physical motivation for such a gen- 
eralization was the possibility of a "gravitational bounce" 



at the late stages of collapse, that could prevent the for- 
mation of a black hole ||. However, a detailed analysis 
by Ori O] showed that shell crossings are inevitable in 
the collapse of weakly (the absolute value of the specific 
charge is less than unity) charged spherical dust shells. 
Inner shells are gravitationally more weakly bound than 
outer shells — due to the larger Lorentz repulsion and 
weaker gravitational potential — and thus collapse slower 
than outer shells, thereby leading to shell crossings before 
the (otherwise) minimum bounce radius is reached. 

One might naively hope that the introduction of a 
positive cosmological constant, A, could prevent the oc- 
currence of shell crossing — thus allowing for singular- 
ity avoidance — by increasing the outward acceleration of 
outer shells with respect to that of inner shells, i.e., by 
making outer shells more weakly bound (since the ef- 
fective 'gravitational potential' contains a repulsive Ar 2 
term that dominates at large radii) relative to inner ones. 

In this paper, we generalize the existing results for 
shell crossing in weakly charged dust collapse, to in- 
clude the contribution of a positive cosmological con- 
stant. Our motivation is fourfold: (i) Recent observations 
of high-redshift Type la supernovae |l^,|lj| and peculiar 
motion of low-redshift galaxies [[l5|, appear to indicate 
that the present radius of the universe is accelerating, 
thus suggesting the existence of a positive cosmological 
constant, A > 0; (ii) such a non- vanishing cosmologi- 
cal constant qualitatively changes the standard asymp- 
totically flat picture of charged dust collapse, allowing, 
in particular, for an altogether different causal struc- 
ture in the static limit: the Reisser-Nordstrom-de Sitter 
(RNdS) spacetime fll6fl ; (iii) the final exterior geometry 
of spherically charged dust configurations is similar to 
that of (more realistic) neutral rotating configurations 
fll7f ; (iv) the inclusion of a positive cosmological con- 
stant could, conceivably, prevent the occurrence of shell 
crossing, thereby allowing — at least in principle — for a 
singularity-free 'bounce' model. 

For definiteness, we shall refer to spherically symmetric 
solutions of the Einstein equations with charged dust and 
a cosmological constant, as generalized Tolman-Bondi 



metrics. 

Contrary to the naive expectation, we find that the in- 
clusion of a positive cosmological constant does not pre- 
vent the occurrence of shell crossing: the latter always 
occurs, irrespective of how large (but finite) the former is, 
for time-symmetric initial data. For non time-symmetric 
initial data, the occurrence of shell crossing is determined 
by the strength of the initial data; for sufficiently steep 
density gradients and non-negative initial velocity gradi- 
ents, a non-zero-measure set of shells will cross at finite 
area radius near the center. 

This paper is organized as follows: Section II derives 
the Tolman-Bondi family of metrics from the Einstein 
equations and discusses shell crossings in such space- 
times. In Sec. Ill, the spherical collapse of weakly 
charged dust with a cosmological constant is studied; the 
Einstein-Maxwell equations are reduced to a coupled set 
of first-order PDE's and an energy-type equation is ob- 
tained. Section IV gives a proof of the inevitability of 
shell crossing, near the center, based on a detailed anal- 
ysis of the shape of the effective potential in the energy 
equation. In Sec. V, we give an alternative proof of the 
occurrence of shell crossing in generalized Tolman-Bondi 
spacetimes, based on the point-particle/dust-shell anal- 
ogy (free surface approximation) . Section VI contains an 
analysis of the uncharged case, with A > 0. Section VII 
concludes with a discussion and summary. 

Natural geometrized units, in which G = c = 1, are 
used throughout. 



II. TOLMAN-BONDI SPACETIMES 

For completeness, we present here a brief description 
of Tolman-Bondi metrics and shell crossings in neutral 
dust collapse in asymptotically flat spacetimes. 

The Tolman-Bondi family of solutions is given by a 
spherically symmetric metric, written here in normal 
Gaussian coordinates {£, r, 6, (f>}: 

ds 2 = -alt 2 + e-™^dr 2 + R 2 (t, r)dSl 2 , (1) 

dn = d9 2 + sin 2 ajJ ~ 2 



(2) 

together with the stress-energy tensor for dust: 

Tab = p(t, r)u a u b = p8 a 5l, (3) 

where u a = Sf is the 4-velocity of a dust element and 
p(t, r) the energy density. 

With the metric ([fl), the independent non- vanishing 
Einstein tensor components are 

G u = i? _2 [-i?e 2 *(2i?'*' + 2R" + R- l R' 2 ) 

-2Rf>R+l+R 2 }, (4) 

Grt = -2R- 1 {R' + R'i>), (5) 

G rr = -R- 2 [e^ 2 * (2RR + R 2 + 1) - R' 2 ] , (6) 



G ee = sin" 2 QG^ = R(M + i?'*'e 2 * 
+i?"e 2 *-i? + *i?-* 2 i?), 



(7) 



where the overdot and prime denote partial differentia- 
tion with respect to t and r, respectively. 
Introducing the auxiliary functions 



2* d'2 



k(t,r) = l-e 2W i? 
m(t,r) = -R(R 2 + k), 

Einstein's equations]] simplify greatly to 

R 2 = 2mR- x - k, 
k = 0, 
rh = 0, 



with the constraint 



to' = 4:TrR 2 R'T ti 



(8) 
(9) 



(10) 

(11) 
(12) 



(13) 



The model is then reduced to a single first-order PDE 
for the area radius R(t, r), given by Eq. (TOT) with 



TO(r)=47r / R 2 {0,f)R'{0,f)p(0,r)df. 



(14) 



If m(r) tends to a constant at spatial infinity, then M = 
lim r ^ +00 m(r) is the ADM mass of the spacetimc. 

The metric (g) with e 2 * = [1 - k(r)]/R /2 , together 
with Eqs. (Eol) and (113), fully determine the Tolman- 
Bondi family of solutions. Included in this family are the 
Schwarzschild metric (m — const.), the Einstein-de Sitter 
universe (R tx rt 2 ' 3 , k = 0), and the closed Friedmann 
universe [R = ra(t), k = r ]. 

The general[j] (k > 0) Tolman-Bondi solution can be 
easily obtained by parametric integration of Eq. (ftm : 



t(rj, r) — io( r ) + fnk 2 (r) + sin?y), 

„/ x 2m „ /"n\ 
R(v,r) = -cos 2 (^, 



(15) 
(16) 



* Since there are only three functions to be determined and 
four equations, only three of these are independent, with the 
remaining one acting as a constraint. We take Eqs. (ho|)-(|l2) 
as our complete set, and Eq. (113) as the constraint equation, 
since it provides a simple relation between the initial data and 
the initial mass profile. 

'The fe > solution corresponds to gravitationally bound 
configurations, and the k — case — discussed below — to 
marginally bound systems. For the unbound case, k < 0, an 
analytical solution has also been obtained in closed form (see 
e.g. |lq] for a systematic treatment of spherical dust space- 
times) , but it is of little interest to the implosion situation we 
are interested in. 



where < 77 < 7r, and to(r) is an arbitrary constant of 
integration to be fixed by the initial data, R(0, r) = v(r), 
via 

t (r) = mi;- 3 tan 3 ( ^M ) [ m (r) + sin m (r)} . (17) 



For time-symmetric initial data, v(r) — to(r) = 0, which 
implies k(r) = 2m/r. 

The radial coordinate r is merely a label for the differ- 
ent shells, and we can therefore fix the radial coordinate 
gauge by equating the initial area radius to the coordi- 
nate radius: 



R(0,r) = r, 
so that Eqs. ( |l5| ) and fll6| ) simplify to 

.3 



t(v,r) = — (77 + sin 77), 
\8m/ 

R(i],r) = r cos 2 — . 



(18) 

(19) 
(20) 



A shell with initial proper area 4irr 2 will thus collapse to 
vanishing area radius in a (comoving) time 



icoii(r) 



(21) 



For inhomogeneous mass distributions (to 7^ const, x r 3 ), 
different shells become singular at different times; in the 
homogeneous case, all the shells collapse to zero area ra- 
dius at the same time . 

The relevant derivatives of the area radius are 

*— 2 (IK (*-£)<" +8 ""MI)- < 22 » 



I? 



tan 



(I)- 



(23) 



A. Shell crossing in Tolman-Bondi spacetimes 

In the context of Tolman-Bondi metrics, shell crossings 
are defined by the loci of events given by 



R' = and R > 0. 



(24) 



At R = R' = 0, a shell focusing singularity is said to 
occur. Unlike shell crossings, this central shell focusing 
singularity does not admit any metric extension through 
it and the spacetime is therefore geodesically incomplete 
(see e.g. jl9|). It has been shown that shell focusing sin- 
gularities can be naked [pOpif l and gravitationally strong 
(finite physical volumes are 'crushed' to zero at the singu- 
larity; see [£2| and references therein) , although massless 



p3| . For the remainder of this paper we shall be con- 
cerned only with shell crossing singularities, as defined 
by©. 

Clearly, the necessary and sufficient condition for shell 
crossings is 



icon < 0, 



(25) 



which is a condition on the initial mass distribution 

to/ 3 
to r 



(26) 



Only sufficiently steep density distributions can provide 
a large enough specific binding energy, k(r), that allows 
outer shells to overlap inner ones. This condition is the 
same as that obtained by requiring R' = 0, with R' given 
byEq. ©. 



1. Marginally bound configurations 

The special case of k(r) = corresponds to a 
marginally bound matter configuration. It is the simplest 
case, and all the relevant expressions can be obtained an- 
alytically in closed form. 

The dynamical equation 



R 2 = 



2to 
~R' 



is trivially integrated to 
R(t , r) = r 



where 



tcW = Vs- 



*c(r) 



2T 3 

9to' 




(27) 



(28) 



(29) 



is the proper time for the complete collapse of a spherical 
shell with initial area radius r. Shell crossings occur when 



R' 



7(r) 



1 



1 



t 



tcY(r) 



0, 



r to 

3 TO 



R' vanishes at 



U 



t c (l+t c ~f). 



(30) 

(31) 

(32) 



Hence, the necessary and sufficient condition for shell 
crossing is 



Oa "^ Or 



(33) 



which is exactly the same as the one for the time- 
symmetric k > case, to'/to > 3/r. [Note, however, 
that the k = case corresponds to non-time-symmetric 
initial data, R{0,r) — ±y/2m/r]. 



III. GENERALIZED TOLMAN-BONDI 
SPACETIMES 

We consider a generalization of the original Tolman- 
Bondi metrics, that includes a charge density distribu- 
tion fJ,(t,r) and a cosmological constant, A. As pointed 
out by Ori |ll| , the existence of a non- vanishing charge 
density implies that the charged matter shells do not fol- 
low geodesic motion, and thus their comoving time t no 
longer equals proper time r, which precludes the use of 
Gaussian normal coordinates (gtt = —1)- We therefore 
consider a general spherically symmetric metric in the 
standard spherical coordinates {£, r, 6, (f>}: 

ds 2 = _ e 2Ht,r) dt 2 + e -2*(t,r) dr 2 + R 2 ^ ^^ ( g 4) 

where R(t, r) is the proper area radius and dft 2 = dO 2 + 
sin 2 9d(jP is the canonical metric of the unit two-sphere. 
The stress-energy tensor is that of a charged dust fluid 
with a cosmological constant, A: 

T a .b = ^-{F ac F b c ~ \g ab F de F de ) + pu a u b - — g ab , (35) 
An 4 87T 

where F ab is the Maxwell tensor, p is the energy density, 
and u a = e~*J° is the four- velocity of a charged dust 
element. F ab is constrained by the Maxwell equations 



\J b F ab = -4irj a , 
Vr F 6cl = 0, 



(36) 
(37) 



where the four-current is j a — pu a , and p is the charge 
density. Spherical symmetry implies that the only non- 
vanishing component of the Maxwell tensor can be chosen 
to be 



Frt-e — 2 , 



(38) 



where Q(t,r) is the total charge inside a shell T, t with 
proper area 47ri? 2 (£,r): 



2n ck pr 



Q(jt t r)= / fdiZ a = / / / 
Js t Jo Jo Jo 



-gj dr'dOdtfi 



Air I p,e'^R 2 df. 

'o 



(39) 



With the metric (p4) the non- vanishing Einstein tensor 
components are: 



G, 



= 2(*+*) 



R 2 

R 

R 2 



[R' 2 + 2R(R" + V'R')] 



(R + 2RV), 



G tr = -(&- R& + tf #) , 



(40) 
(41) 



Ct ,■ 



,-2(*+$) 

_ R 2 



2RR + R(R - 2i?$) 



~(# + 2R$% 



(42) 



= -R{( 



-2* 



-R + i?($ + %) + R(V - * 2 - *$) 



+e 2 * [R" + #'($' + *') + i?($" + $' 2 + $'*')] }> (43) 
G^ = sin 2 6Gee. (44) 



The non-vanishing components of the stress-energy 
tensor arc 



e 2 * / ' G 2 

Ttt = — ( ^7 + A + 87T/7 



8tt Vi? 4 

n -2* 



T = 



87T 



i? 4 

i? 2 



Q 2 
I? 4 



Tee = sin flT^ = — - -A + — 
on \ 

Introducing the auxiliary functions 
k(t,r) = l-e 2 *R' 2 , 



the Einstein-Maxwell equations can be recast as 

A 



-2$ n2 



R z = -Q Z R- Z + 2mR- L -k+ -i? 2 , 

k = 2 R \(k-l)&, 

rh = 0, 

Q = o, 

rri = 4irR 2 R'p + QQ'R- 1 
Q' = AirR 2 e-*p. 



R 2 R' 



(45) 
(46) 
(47) 

(48) 
(49) 

(50) 

(51) 

(52) 
(53) 

(54) 

(55) 



Equation @ is just Eq. @; Eqs. (P)- ^) fo llow di- 
rectly from the Einstein equations, and Eq. (|55j) follows 
from Eq. (B9[). The equation for k' is algebraically com- 
plicated and it proves more convenient to use an equiv- 
alent equation given by the local conservation of energy- 
momentum, V a T fc a = (where we took b = r): 



QQ' = AirpR 4 ®' 



(56) 



This equation expresses the balance between Lorentz and 
gravitational 'forces' in the comoving frame. The other 
equation (b = t) from the local conservation of energy- 
momentum gives the conservation of specific charge: 



d f . I £ I = e = 0. 



(57) 



Equations ([48|), (|50|), (pTj) , and ( |5^ ) form a complete set; 
the remaining independent equations are constraints. We 
have thus reduced the Einstein-Maxwell system to a set 
of coupled first-order PDE's. It is not the purpose of 



the present paper to obtain an explicit general solution 
in closed form — this has been done, in a different co- 
ordinate system, by Ori O). As we shall see shortly, 
first integrals of the equations of motion suffice to prove 
the inevitability of shell crossings in generalized Tolman- 
Bondi model spj 

Using Eqs. (|§), (§|), and (f|), Eq. (|l|) integrates to 

1 2 



k(R,r) = 1 



R 



W(r) 



(58) 



where W(r) is a free function to be fixed by the initial 
data, via Eq. (|30|) evaluated at t = 0. Comparison with 
the neutral case, leads to the interpretation of W as a 
measure of the binding energy of the system per unit 
mass. For A = 0, the configuration is gravitationally 
unbound, marginally bound, or bound, when W > 1, 
W = 1, or W > 1, respectively. 

For time-symmetric initial data, -R(0, r) = 0, with the 
scaling i?(0, r) = r, we have then 



W(r) = e— 

r 



2m 
1 + ^ 



2 



3 r 



(59) 



where the '+' sign was fixed by consistency with the 
asymptotic A = Q = limit. 

Regularity at the origin requires p0| : 



lim — =0, 

r^0+ r 


(60) 


to 
lim — = 0, 

r^0+ r 


(61) 


lim |e| = |e |e(0,l), 

r->0+ 


(62) 



which implies 

W(0) = W = 1. (63) 

Summarizing, the only non-trivial dynamical equation 



e-^R z + U(R, r) = 0, 



(64) 



with the effective potential 




U(R, r) = a + aiR^ 1 + a 2 R~ 2 + a 3 R 2 , 


(65) 


oo(r) = 1 - l^ 2 , 


(66) 


oi (r) = 2m(e— W-lj , 


(67) 


a 2 (r) = Q 2 (l-e 2 ), 


(68) 


A 
a3 = "3' 


(69) 



"'"This is essentially due to the fact that shell crossings are 
dynamical processes, whence the knowledge of the first 'time' 
derivative of the spheres' radii suffices to analyze their relative 
motion. 



IV. EFFECTIVE POTENTIAL ANALYSIS 

In this section, we prove two Propositions that show 
that shell crossing is inevitable in (weakly charged) gen- 
eralized Tolman-Bondi spacetimes. 

Proposition 1. For weakly charged initial data, there 
exists r* > 0, such that the effective potential U(R,r) 
has, for r e (0,r*), three distinct real positive zeros, 
Ri < Ri < i?3, satisfying (see Fig. 1) 



C/(i?i <R<R 2 ,r) <0, 
U(R 2 <R<R 3 ,r) >0. 



(70) 

(71) 



Sketch of the proof. We prove the above proposition 
by showing that for < r < r* <C 1, U(R,r) has a lo- 
cal negative minimum at R m > and a local positive 
maximum at Rm > Rm- Since \ira^^ Q + U = +oo and 
limij_ > + 00 U = — oo, it follows that U(R, r) has three dis- 
tinct real positive zeros, obeying conditions J7Q)-(71). 

Proof. Let us rewrite 



U(R,r) =a 3 R- 2 H{R,r), 

H(R, r) = R 4 + c 2 R 2 + c x R + c , 

co^-jQHl-e 2 ), 



C\ 



-m eW— - 1 
A V to 



c 2 ^-^(l-W 2 ). 
The relevant limits are now 



lim H = c < 0, 

lim H = +oo. 



(72) 

(73) 

(74) 
(75) 
(76) 

(77) 
(78) 



Although the local extrema of U are different from those 
of H, all (non-zero) roots of H are also roots of U, and we 
can thus study the qualitative behavior of U by analyzing 
H and its first partial derivative with respect to R: 



IdH 
TdR 



R A 



1 1 

2 C2i?+ 2 Cl - 



Now, djiH has three distinct real zeros provided 



W 2 -l 



<-A 



■ Q 



6to I eW 1 

m 



-I 2 



(79) 



(80) 



Since the right-hand-side of the inequality is manifestly 
negative, this requires W < 1. Taylor expanding Q(r) 
and m{r) near the origin, using Eqs. (54)-(J56[), and Eq. 
( pq ) evaluated at t = 0, with the scaling i?(0, r) = r, we 
have 



W(r) = 1 



4~ 



Po(l 



0(r 4 



(81) 



for r£ (0, £1), where 5i « 1. Condition ( J80| ) is trivially 
satisfied for < r < £i. 

We will now show that of the three real roots of 
8rH, two are positive, say R\ and R 2 , such that 
H(Ri)H{R 2 ) < 0, thereby guaranteeing that H has a 
local positive maximum at R\ £ (i?i,i^j) and a local 
negative minimum at R 2 € (R2,Rs)- 

The roots i?i (i = 1, 2, 3) of the polynomial 9fli? satisfy 



$^ = 0, 






(w 2 



2A 

2A V to 



and can be explicitly given as 

'e 

'e 2tt 



_Ri = 2-^/xcos 
R 2 = 2^008 



where 



R 3 = 2Vxcos ( - + — 



_l-W 2 

X ~ 2A ' 

/A eQW - m 
cos 6 = S\ — - 



(82) 
(83) 

(84) 

(85) 
(86) 
(87) 



v 2 {1-W 2 ) 3 / 2 ' (89) 

Since cos is anti-periodic in 7r and the trigonometric ar- 
guments in Eqs. (85)-(p6|) are separated by a 2-7r/3 phase, 
it sufhces to require that two of the roots are positive, the 
third one being then necessarily negative. We therefore 
impose 7r/2 < < 37r/2 (thereby guaranteeing that R\ 
and R 2 are positive), which translates to 



3m(l-e^|)< V f|(l-W a )». 

Near the origin, for < r < ^ < 1, this becomes 



(90) 



47F f, 2N A 

T Mi-^) + - 



r*+O(r 3 )>0, (91) 



which is trivially satisfied; thus R\ and R 2 are positive, 
for r e (0,£ 2 ). 

Now, it remains to show that H(Ri)H(R 2 ) < 0. From 
Eq. ( f73| ) together with properties (p2j)-(p^), and setting 
i?3 = —x, we obtain 

1 /C2N2 

"4V4 



HiRJHfa) = f(x) 



2 -2 



1 . /Ci 

i C ClC2 + (j 



+^ + |(j) - 



-CqCiX 



(92) 



Using Eqs. ([74[) - (f76|) , it is straightforward to check that 
/ is a monotonically increasing function of x. From Eq. 
( p7| ) we have x < \/—c 2 /6 = %. Thus, it suffices to 
show that /(%) < 0: 




C\ 



G 

-C2 



<0. 



(93) 



Near the origin, for r £ (0, £3), with £3 <C 1, we have 



G 



-OLi 



(^)V + 0(r«), 



where 



6 /4tt A\ . 

ai = AlT Po+ 3' (1 " eo) ' 



«2 



3 

A 



A-n , 9 , A 
yPo(l-£o)+ 3 



(94) 

(95) 
(96) 



Condition ( |93| ) is manifestly satisfied; hence 
H{R 1 )H{R 2 ) < 0, for < r < &. 

Now, let us take r» = min {£1 , £2 , 6s}_, such that for 
r e (0, r») conditions (p0|), (|9C|), and ( |93| ) are mutually 
satisfied. 

Summarizing, we have shown that there exists r» > 0, 
such that for r £ (0, r»), d^if has two distinct real posi- 
tive zeros, corresponding to the local extrema of H (R, r). 
We have further shown that the product of the local ex- 
trema of H is negative, which implies that H has three 
distinct real positive roots, since lmiR^o+ H = cq < 
and lim^^ +oc H = +00. By construction, these are the 
same zeros of U(R, r). Since U(R, r) — a 3 (r)R~ 2 H(R, r) 
and 03 < 0, it follows that U(R, r) is negative for 
R £ {R\,Ri) and positive for R £ (R 2 ,R 3 ), as desired. 
■ 

Proposition 2. Consider a general solution of Eq. 
(pp|), with parameters {Q, to, A} satisfying Proposition 1, 
and the spacetime metric determined by that solution. In 
such a spacetime, shell crossing always occurs from reg- 
ular time-symmetric initial data, for a non-zero-measure 
set of spherical weakly charged matter shells. 

Proof. For r > 0, \im R ^ a + U(R,r) = +00 and 
linif{^ +00 U(R,r) = —00. Since U is C°° with respect 
to R, this implies — from Eq. (pGJ ) — that there is a time 
t m > such that 



R(t m ,r) = 0. 



(97) 



Let us then define Rm(r) = R(t m ,r). 

Taking the total derivative of Eq. ( pCf ) with respect to 
r, we obtain 

. . . r)TT 

2e~ 2 *RR' ~ 2$'e- 2 * J R 2 + ^-R 1 + U' = 0. (98) 

ait 



For time-symmetric initial data, R(0, r) = 0, and thus, 
from Eqs. ©-© we have 



R'(0,r) 



-< T '^<(m 



R=r 



R'(t m ,r)=-U'(R m ,r)/(^j 



(99) 
(100) 



We now examine the signs of the numerators in the equa- 
tions above. From Eq. (Q) we have 

U'(r,r)=a' + a' 1 r- 1 +a' 2 r- 2 , (101) 

U'{R m ,r) = a' + a'iR- 1 + a' 2 R^ 2 . (102) 

Setting R(0,r) = r = R 2 and R m (r) = R\, it follows 
from Proposition 1 that R m (r) < r, and we thus can set 
R m = P r r for any given shell r, where /3 r < 1 depends 
on r. Near the origin, for < r < r* <C 1, we have then 



U'(r,r) = 
U'(R m ,r) = 



yPo(l-e2)+A(!-e8 



+0(r 3 ), 
6 

Jr' 

A 



47T 

4 I yPo(l 



, v l-e^+3/3 r )]r + 0(r 3 ). 



(103) 



(104) 



Both numerators are manifestly negative, for r £ (0, r^.). 
Now, from Proposition 1, we have (dU/dR)u =r > 
and (dU/dR) R=Rm < 0. Since 17(0, r) = E7(t m ,r) = 0, it 
follows from Eqs. (|9S|)-(|T00|) that 

R'(0,r) 



R'(t m ,r) 



<0. 



(105) 



Hence, by continuity, there is a time i sc £ (0, i m ), such 
that 



#(*«,, r) = 0, Vre (0,min{r*,r*}). 
This completes the proof. ■ 



(106) 




FIG. 1. Effective potential U(R,r) and auxiliary poly- 
nomial H(R,r) oc R 2 U(R,r), for r £ (0,r*). With- 
out shell crossing, a spherical shell with initial area radius 
R(0, r) — r = R2 would collapse towards decreasing values 
of R until reaching the minimum bounce radius H m = Ri at 
some later time i m . However, shell crossing is inevitable at 
i sc £ (0,i m ). 



V. FREE SURFACE APPROACH ANALYSIS 

This Section gives an alternative proof of the in- 
evitability of shell crossing in spherical weakly charged 
dust collapse with a cosmological constant. It relies on an 
immediate consequence of the generalized Birkhoff theo- 
rem p5fl (which can be extended to include a cosmolog- 
ical constant p6[), the so-called free surface approach: 
as long as shells do not cross, their motion will depend 
solely on their total interior charge and mass, and thus 
each shell will move independently (albeit not geodesi- 
cally) like a test particle in a RNdS background, with 
mass and charge equal to the total mass and charge of 
that shell (see e.g. ||). 



A. Particle motion in a Reisser-Nordstrom-de Sitter 
geometry 

In Schwarzschild coordinates {t, r, 8, </>}, the RNdS 
metric reads 



ds 2 = -Adt 2 + A^dr 2 - r 2 dn 2 , 

A r =1 + -j- - -r , 

r r z 3 



(107) 
(108) 



where dfl 2 = d9 2 + sin 2 6dtp 2 is the canonical metric of 
the unit two-sphere. 

The action for a particle of mass rh and charge q is 

S[x a , e] = JdX (~e- l x a x h g ab - ^fh 2 e - qx a A^\ , 

(109) 

where A is an affine parameter along the particle's world- 
line, A a — (Q/r)<5£ is the electromagnetic four-potential, 
and e = e(A) is the einbein, an independent function that 
generalizes the action for a point particle J27J ]. 

It is straighforward to show (cf. Appendix) that if 



(£tA) a = ?A a , b + M% = 0, 



(110) 



for a Killing vector field £ a , then S is invariant, to first- 
order in £ a , under the transformation x a — * x a + at; a , 
where a £ R. Denoting by 7r a = rhu a — qA a the general- 
ized momentum of the particle, it then follows that the 
corresponding Noether charge, £ = £ a 7r a , is a constant of 
motion: 



d£ 

rfA 



(Tr a VbC+eV b ir a )u b 



,a n b 



V 6 £ a + q(C^A a )u a = 0. 



(Ill) 



One of the trivial Killing vectors of the RNdS metric 



is £(t) = d t - Thus, with the metric (107), we have 



E = t?(mut - qA t ) = mA— + q—, 

dr r 



(112) 



where r is the particle's proper time. Since the geometry 
is not asymptotically flat, the constant of motion E can- 
not be interpreted as the particle's total energy measured 
at spatial infinity; it should be regarded as an effective 
energy, that coincides with the particle's total energy at 
infinity in the A — > limit. 



From Eq. (107) we have, for a timelike radial worldline 

2 / , \ 2 



-I 



1 



*-•£ 



which, together with Eq. (112), gives 



+ V(r) = 0, 



V(r) = (1 - W 2 ) - 
Mfi -2^ 2 A 



I -We 



Q_ 

' M 



2M 

r 



(113) 



(114) 



(115) 



where W = E/m is the particle's specific effective energy, 
and e = q/fh its specific charge. 



B. Shell crossing analysis 



Equation (114), which governs the radial motion of 
test particles with specific charge e = q/fh in a RNdS 
geometry with fixed total charge Q and mass M, will co- 
incide with Eq. (p4), which governs the radial motion of 
charged dust shells with total interior charge Q(r) and 
mass m(r), provided we make the following formal iden- 
tifications [and set $ = in Eq. ( p4| ) , thereby identifying 
proper and comoving times]: 



r = R, 

W = W, 

M = m, 

Q = Q, 

- Q 



(116) 
(117) 
(118) 
(119) 

(120) 



geometry, whereas the former travels on a geometry that 
is determined by the interior mass and charge of each 
shell. 

Let us now consider two infinitesimally close shells, 
with wordlines ro(r) and T\(t) = ro + £ (where < £ <C 



1), which are solutions of Eq. (114) with V(r) < 0. As 
shown in Fig. 2, these two neighboring shells will cross 
provided 



dr\ 



dr\ 
dT~) 



(122) 



ro+i 



From Eq. (114), this translates into V(r + £) < V(tq). 
Expanding V(r Q +£) = V(r ) + V'(r )£ + C(£ 2 ), to first 
order in £ we have then V'(ro)£ < 0. Since £ > 0, the 
condition for shell crossing is simply (dropping the '0' 
index) 



V'{r) < 0. 



(123) 



Evaluating Eq. (|12l|) near the origin, for < r < r c <C 1, 
this condition becomes 



2e 2 A(l + J-] r + O(r 3 )>0, 



1+ 

8-KpoJ 



(124) 



which is clearly satisfied for r £ (0, r c ). Shell crossing is 
inevitable near the center. 




(dx/dr) 



The potential V(r) reads then 

V(r) = (l-W 2 )-(l-We®) — 
\ m I r 



1 



Q 2 \Q 2 A 



m- 



-r 2 . 
2 3 



(121) 



It governs the radial motion of a charged dust shell with 
proper area 47rr 2 and specific charge e{r) = Q(r)/m(r), 
in a RNdS background metric with total mass m(r) and 
charge Q(r). In the free surface approach, the only — but 
crucial — difference between shell and particle motion in 
a RNdS background, is that the latter travels on a fixed 



FIG. 2. Shell crossing for spherical weakly charged dust 
collapse inr-r coordinates. For ro G (0, r c ), any two neigh- 
boring shells, ro and ro + £, will cross each other at r = r sc , 
when V'(ro) < 0. 



VI. NEUTRAL CASE 

In this Section, we analyze separately the occurrence 
of shell crossings for neutral dust spheres in the presence 
of a positive cosmological constant, which is a simpler 
and more realistic model for macroscopic astrophysical 
objects. 



Since there is no elecrostatic repulsion acting on the 
test particles, they move geodesically in a Schwarzschild- 
de Sitter geometry |29]. We can therefore set g tt = — 1, 
thereby identifying proper and comoving times. The rel- 
evant dynamical equation is then 



U(R,r) 



1 



R 2 + U{R,r) = 0, 

9 2m A .-, 

W -It- 3 R - 



(125) 
(126) 



We consider here the particular case of W = 1, cor- 
responding to a gravitationally unbound configuration, 
since it allows for an analytical solution in closed form: 



""•'■' = ('x 



9(t,r) 



sinh 2/3 6(t,r), 



[*c(r)-t] 



(127) 



(128) 



where t c (r) is an arbitrary function to be fixed by the 
initial velocity profile via 



R(0,r) = - 



coth 



7 3At c 



The relevant derivatives of the area radius are 
R'(t,r) =R(^- + J^t' c cothe\, 



R(t,r) = - 



R coth 9. 



(129) 



(130) 



(131) 



Since 9 coth 9 > 0, the necessary and sufficient condition 
for shell crossings is 



t <0. 



(132) 



A trivial example is t c — 1/r. For small r, 3m/m' ~ r, 
and coth 9 ~ 1/9. Solving for R' = near the origin 
yields then 



t sc (r) 



1 
3^ 



<t c 



(133) 



thus confirming the occurrence of shell crossing near the 
center. 

Since the initial velocity profile [cf. Eq. (129)] is a 



monotonically increasing function of t c (r), the condition 
t' c < is simply the requirement that outer shells have 
a sufficiently larger inward initial radial velocity than in- 
ner ones, to overcome the (comparatively larger, for outer 
shells) repulsive effect of A, thereby leading to shell cross- 
ing before complete collapse. 

We note that since W — 1 corresponds to unbound 
systems, the above criterion for shell crossing is likely 
to be relaxed for bound configurations. For W < 1, for 
each shell r, there is a critical value A c (£) — which is the 
root of Eq. (126) for fixed r — such that the system is 
gravitationally bound for A < A c , and unbound for A > 
A,. 



VII. DISCUSSION AND CONCLUSIONS 

We have shown that the inclusion of a positive cosmo- 
logical constant in spherical charged dust collapse does 
not prevent the occurrence of shell crossing near the cen- 
ter. Heuristically, this can be explained by the fact that 
A has a long range effect, whereas the relevant physics 
for shell crossing occurs near the center, where Lorentz 
and gravitational 'forces' compete and the A repulsion is 
negligible. 

Although the free surface approach yields a rather sim- 
ple method for proving the inevitability of shell crossing 
near the center, being a purely kinematical analysis it 
cannot relate the initial data to the motion of the shells. 
Reducing Einstein's equations to first integrals of motion 
enabled a more physical analysis, where the dynamics of 
collapse is determined by the choice of initial data. 

Proposition 1 showed that in a small neighborhood 
near the center of the matter distribution, the effective 
potential is such that a shell with initial area radius 
r = i? 2 will collapse towards smaller values of R and, 
provided there is no shell crossing, reach the minimum 
bounce radius R\. Proposition 2 showed that for a such 
a shell [R(0, r) = r = R2] shell crossing will inevitably 
occur before the minimum bounce radius is reached. 

Proposition 2 used the assumption of time-symmetric 
initial data. This ansatz is by no means essential. For 
non-time-symmetric implosion situations, two possibili- 
ties arise: (i) the initial data is sufficiently strong, such 
that all shells have R(Q, r) < 0, or (ii) a non-zero-measure 
set of shells has R(0, r) > 0, and therefore initially starts 
to expand towards increasing values of R, until each shell 
reaches a maximum area radius -R m ax(?") and then col- 
lapses back through its initial radius. Case (ii) reduces 
to the time-symmetric situation when i?(£, r) = i? max 
and is thus covered by Proposition 2. For case (i), the 
occurrence — or lack thereof — of shell crossing is deter- 
mined by the initial velocity profile, v(r) = i?(0,r) for a 
given mass m(r) and charge Q(r) [hence specific charge, 
e(r)] distribution. Clearly, if v < and v' = then shell 
crossing will occur, since it does when v = v' = 0: all the 
shells are differentially accelerated in the same manner, 
irrespective of their initial velocity profile. It then follows 
that shell crossing will occur for any v(r) < 0, provided 
v 1 < 0. 

For the neutral case, shell crossing occurs — even 
for gravitationally unbound matter configurations — 
provided the initial velocity profile is sufficiently steep, 
irrespective of how large (but finite) A may be. Unlike 
Lorentz 'forces', which become more noticeable as col- 
lapse proceeds and the area radius of the shells decreases, 
the A repulsion becomes increasingly unimportant at late 
times (i.e., in the strong-field region, at small radii), 
whence the criterion for shell crossing becomes analogous 
to that for uncharged dust collapse in an asymptotically 



flat spacetime. 
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APPENDIX 

Consider the following action: 

S[x a ,e] = / dX ( -e~ 1 x a x h g a b - -m 2 e - qx a A a 

(134) 

where A is an afflne parameter along the particle's world- 
line, A a is the electromagnetic four-potential, and e = 
e(A) is the einbein. Assume that there is a Killing vector 
field £ a , such that (£^A) a — 0. Under the transformation 
x a — > x a — x a + er£ a , where a € R, to first-order in a we 
have 



(135) 
(136) 

(137) 



S= J d\[-e-\i a + aZ%x c ){± b + ae, c x c ){g ab + ag ab , c C) 



x a -» x = x a + a£ a c x c , 

9ab ~ > 9~ab = 9ab + CT9ab,c^ C , 
A a ► A a = A a + uA a ^ c t, ■ 

The action S[x a ,e] becomes then 
r l 



1 



TO 2 e _ 9(± a + ^ x cp a + aAa ^c 



= S + ^JdX [e- 1 x a x\C i g) a b - 2(£ ? A) c f ] 



S. 



(138) 
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